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The magnetic field dependence of the electronic properties of real single vertical quantum dots 
in circular and rectangular mesas is investigated within a full three-dimensional multiscale self- 
consistent approach without any d priori assumptions about the shape and strength of the con- 
finement potential. The calculated zero field electron addition energies are in good agreement with 
available experimental data for both mesa geometries. Charging diagrams in a magnetic field for 
number of electrons up to five are also computed. Consistent with the experimental data, we found 
that the charging curves for the rectangular mesa dot in a magnetic field are flatter and exhibit less 
features than for a circular mesa dot. Evolution of the singlet-triplet energy separation in the two 
electron system for both dot geometries in magnetic field was also investigated. In the limit of large 
field, beyond the singlet-triplet transition, the singlet-triplet energy difference continues to become 
more negative in a circular mesa dot without any saturation within the range of considered magnetic 
fields whilst it is predicted to asymptotically approach zero for the rectangular mesa dot. This dif- 
ferent behavior is attributed to the symmetry "breaking" that occurs in the singlet wave- functions 
in the rectangular mesa dot but not in the circular one. 

PACS numbers: 73.21.-b 



I. INTRODUCTION 



In recent years the electronic properties of semicon- 
ductor quantum dots (QDs) have attracted considerable 
attention because these zero-dimensional objects exhibit 
properties characteristic of atomic few-electron systems. 
QDs are also regarded as promising candidates for future 
device applications [lj, not least of which is the possi- 
bility of realizing a scalable quantum computer 0|. A 
considerable amount of theoretical and experimental re- 
search has been devoted to understanding the electronic 
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properties of cylindrical mesa structures within which the 
QD has a shape of a thin disc. These disc-shaped QDs 
are usually formed from double barrier structures made 
in to small mesas and surrounded by a gate electrode. 
The number of electrons is fully tunable starting from 
zero 0. In such systems three-dimensional quantization 
coupled with rotational symmetry leads to the existancc 
of a "shell structure" in the electron addition energy, i.e. 
for certain electron numbers the electron addition ener- 
gies show pronounced maxima. Spin effects also lead to 
the appearance of smaller secondary peaks in the mid- 
shell regions, in accordance with Hund's first rule 

In contrast to real atoms, electrons in vertical QDs 
occupy a sizable region of space of the order of 100 nm. 
This allows utilization of a relatively small and readily ac- 
cessible magnetic field to probe their electronic structure 
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and also to induce transitions between the various elec- 
tronic states as the number of particles and the strength 
of field are varied. However, perfect circular symmetry 
of the QD is difficult to achieve and maintain in real scal- 
able electronic circuits. Ideal atomistic features are also 
very sensitive to any kind of disruption in the system [f| . 
This dictates the necessity of understanding the response 
of the electronic structure to various geometric deforma- 
tions. 

Compared with circular QDs which have been studied 
by a variety of methods Q, QDs in rectangular mesas 
(hence forth called rectangular QDs for brevity) have at- 
tracted considerably less attention both experimentally 
and theoretically. Invariably, most calculations assume a 
pure two-dimensional confinement potential. Frequently 
El H- El fioj - a model two-dimensional anisotropic har- 
monic oscillator potential (l/2)m*(oj 2 .x 2 + uiyy 2 ) charac- 
teristic of an elliptic QD with confinement strengths lo x 
and u)y in the x and y directions respectively, is used 
in theoretical calculations for this purpose. Simulations 
conducted with such a potential 6] showed that it is not 
straightforward to match the experimental geometric as- 
pect ratio (3 = L/S, where L(S) is the length of the 
longest (shortest) side of the top mesa contact, to the 
potential deformation ratio 6 = Lj y /uj x . Breakdown of 
Hund's first rule was also found in QDs with a deforma- 
tion ratio S as small as 1.2 Evolution of the excita- 
tion spectra in circular and elliptic QDs in a magnetic 
field up to 4 T was also studied 0|- Whilst the pre- 
vious works gave some understanding and clarification 
of the shell structure evolution in response to potential 
deformations, the detailed analysis and comparison with 
experimental data is still lacking. 

In this work we try to remedy this shortcoming and 
present calculations of the zero field addition energy spec- 
tra and charging diagrams in an applied magnetic field for 
real circular and rectangular vertical QD mesa structures. 
We perform a full scale three-dimensional (3D) analysis 
of the QDs without any d priori assumptions about the 
shape or strength of the confinement potential. We uti- 
lize density-function theory (DFT) to describe the many 
body interactions amongst the electrons in the QD, and 
employ a semi-classical Thomas-Fermi approach to treat 
charge evolution outside the QD region. In this multi- 
scale approach, the confinement potential is obtained di- 
rectly from the self-consistent solution of the Poisson and 
Schrodinger (Kohn-Sham) equations with corresponding 
device boundary conditions. Our calculations are based 
on real experimental structures for which well resolved 
addition energy spectra and magnetic field characteris- 
tics have already been measured |3,EI 0- Previously this 
full scale method was applied to a variety of problems 
such as, e.g., the simulation of four-gated QD mesas with 
a square geometry ^ which showed that elliptical defor- 
mations have little effect on the addition energy spectra 
but affected the charge configuration greatly. Prelimi- 
nary calculations for rectangular QDs also showed that 
the shell structure is completely destroyed for f3 > 1.2 



and a Wigner molecule in the two-electron system is 
formed for (3 > 3 [T^ . 

The paper is organized in the following way. The QD 
structures utilized in this work are briefly described in 
the next Section. A detailed description of our model is 
given in Section IlIII In Section llVI we first compare ad- 
dition energy spectra and shell structures of the circular 
and rectangular QDs in the absence of a magnetic field 
( Section flV Aft . Charging processes and evolution of spin 
states under application of magnetic field are then con- 
sidered in Section IIV Bl Comprarison with experimen- 
tally measured charging diagrams in a magnetic field is 
also performed in the same section. The particularly in- 
teresting and important case of a two electron system is 
addressed in Section IIV CI Finally, Section contains 
concluding remarks. 



II. DEVICE STRUCTURES 

Our QD structures have a layered profile in the verti- 
cal z-direction Q which are the same for both the cir- 
cular and rectangular mesa structures. Two layers of 
Alo.22Gao.78As are used as potential barriers for electrons 
confined in a quantum well made of Ino.05Gao.95As. Con- 
finement of electrons in the lateral xy-plane is achieved 
by applying a bias Vg to the gate on the side of the mesa. 
Full details on the geometry, layer composition of these 
structure, and the measurement techniques can be found 
elsewhere HH0 E^. 

For the simulations described below, we represent the 
mesa and substrate as a cylinder or a rectangle with total 
height of 350 nm. The side gate width is taken to be 
equal to 150 nm, while the center of the quantum well is 
located 116.5 nm from the top of the side gate. The lower 
100 nm of the simulation cylinder is assumed to be buried 
in the substrate. The diameter of a circular mesa is equal 
to 0.5 urn in accordance with experimentally measured 
mesa [3J. For the rectangular QD mesas, L and S are 
taken to be 0.55 ^m and 0.4 ^m {(3 = 1.375), and 0.65 
/im and 0.45 fim (/3 = 1.44). These mesas correspond 
respectively to structures identified as rectangular QD X 
and rectangular QD Y studied earlier at low magnetic 
fields @. Accordingly we refer to them as rectangular 
QD X and rectangular QD Y in this paper. 



III. MODEL AND SIMULATION DETAILS 

Within the spin-density functional theory, which we 
use to describe electrons confined in the QD region, the 
charge density is calculated after solving separate single- 
particle Kohn-Sham equations for electrons with spin up 
(t) and down (j): 

fl-TU)^TU)( r ) =£ Ta)^T(i)( r ) j (x) 
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where the Hamiltonian j s given as follows: 

fl-T(i) = f + 0( r ) + A£; c ( r ) + V UD ( r ). (2) 

In this equation T is the kinetic energy operator which 
in the presence of an external magnetic field B reads as: 



f=l(-ihV+-A) — V- 
2 V c / m*(r) 



-jSV + - A 

c 



(3) 



where m*(r) is the position dependent electron effective 
mass, and A = (1/2) (By, —Bx) is the corresponding vec- 
tor potential. We assume here that the magnetic field is 
oriented along the z-direction and neglect Zeeman split- 
ting for clarity. 

In Eq. J2|, AE c (r) stands for the conduction band off- 
set between the different materials. The values of AE c (r) 
are fixed at 180 and -40 meV for Alo.22Gao.7sAs/GaAs 
and Ino.05Gao.95As/GaAs respectively. Note that for the 
real structures with non-zero doping in the GaAs contact 
regions, the value of AE c (r) for Ino.05Gao.95As/GaAs is 
difficult to determine precisely. 

The exchange-correlation potential v xc (r) is com- 
puted within the local spin density approximation 
(LSDA) [13| • Note that this functional does not explicitly 
depend on magnetic field. Comparison of DFT results 
with calculations using current-spin density functional 
theory 14] for two-dimensional systems showed that this 
approximation is reliable over a wide range of magnetic 
field, although at higher fields, effects of paramagnetic 
currents in v xc (r) become more important. 

The potential <f>(r) = (j> ext (r) + </w(r) + (for (r) in Eq. 
@ is the sum of the external potential <fi ex t(r) due to 
the applied voltage, screening potential 4>i on (^) arising 
from the ionized impurities in the structure, and Hartree 
potential 4>h(y) accounting for the repulsive electron- 
electron interactions. It is obtained from the solution 
of the Poisson equation: 



Ve(r)V<£(r) = 4rrp(r), 



(4) 



where e(r) is the dielectric constant of the medium, and 
p(r) is the charge density which inside the QD region is 
equal to 




V> T (r)| 



E 

occup 



t(;^(r) I 



(5) 



with the summations spanning occupied states for elec- 
trons with spin up and down (the number of those 
states is, in general, different). Outside the region, 
charge distribution is determined from electron n(r) and 
hole p(r) densities calculated within the semi-classical 
Thomas-Fermi approximation [TEl Il6| screened by the 
ionized donors and acceptors with concentrations A^(r) 
and JV7(r) (which were assigned their nominal values 



p(r) = e [N+(r) - N A {r)+p{v) - n(r)} 



(G) 



Since vertical QDs are usually much smaller than the 
physical dimensions of the device, the Kohn-Sham wave- 
functions actually vanish long before reaching the device 
boundaries. This allows us to embed a local region in the 
global mesh for solving the Kohn-Sham equations. This 
local region is chosen to be large enough to ensure van- 
ishing wavefunctions on its boundaries. For the Poisson 
equation zero normal electric field on the lower part 
of the simulated structure buried in the substrate and 
on the top contact region are used as a boundary con- 
dition. On other surfaces, not covered by the side gate, 
the potential 4> is set equal to the Schottky barrier value 
Vs = 0.9 eV. Boundary values of the potential on the 
side gate are equal to the Schottky barrier value modi- 
fied by the applied gate bias, Vs —Vg- The system of 
Kohn-Sham and Poisson equations (@J is solved iter- 
atively until a self-consistent solution for the Kohn-Sham 
orbitals ^(^(r) and eigenvalues is obtained. 

The calculations are performed on a parallel platform 
by means of the finite element method (FEM) with tri- 
linear polynomials on a variable size grid [TtI ] . Previo usly 
FEM was applied in electronic structure calculations [18| . 
The advantages of FEM utilization are the ability to sys- 
tematically improve the accuracy by expanding the basis 
set and its perfect suitability for parallclization. 

In our approach Poisson's equation is solved by means 
of the damped Newton-Raphson method while the gen- 
eralized eigenvalue problem obtained after discretization 
of the Kohn-Sham equation is tackled by means of a sub- 
space iteration method based on a simple Rayleigh-Ritz 
analysis |17| . The small number of required eigenpairs 
(in the present work, less than 10) made this approach 
sufficient. A parallel conjugate-gradient method precon- 
ditioned with block Jacobi with an incomplete LU fac- 
torization on the blocks available in PETSc (Portable, 
Extensible Toolkit for Scientific Computation) [l!j is uti- 
lized for solving the resulting matrix equations. In the 
presence of a magnetic field, the matrix obtained from 
the Kohn-Sham equation is Hermitian and the hermitian- 
conjugate method 19] with the same preconditioner as 
above had to be used in solving the eigenvalue problem. 
Compared to the ordinary conjugate gradient method 
with a Jacobi preconditioner, we found that this ap- 
proach gives rise to at least an order of magnitude in- 
crease in performance especially when dealing with Her- 
mitian matrices |2(J. 

After the eigenvalues are determined, the electron 
charging diagram can be calculated. The experimentally 
measured charging diagram is obtained by single electron 
transport spectroscopy HH0, when the QD is in the 
regime of Coulomb blockade with only a very small bias 
applied across the QD in the vertical z-direction. Con- 
sequently, the difference between the chemical potentials 
of the source and drain contacts is negligible and both 
of them can be set to zero without loss of generality. A 
charging event occurs when the chemical potential of the 
QD becomes equal to that of the source (or drain) con- 
tact. Then we can determine the stable electronic config- 
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uration by directly comparing the total energy E(N) of 
the iV-electron system with the system containing N — 1 
electrons, the difference of which gives the required value 
of the chemical potential n(N). Furthermore, we use 
Slater's rule [2lll22| (transition state technique) in order 
to calculate n(N): 

H(N) = E(N) - E(N -1) w e(N -1/2), (7) 

where e(N — 1/2) is the eigenvalue of the state with half 
occupancy in the system with N — 1/2 electrons (the 
transition state). If fi(N) < 0, then the iV-electron con- 
figuration is stable, otherwise the number of electrons in 
the QD is TV - 1. The voltage V G (N) at which fj,(N) = 
gives the charging voltage for the 7V-th electron (or equiv- 
alently it defines the boundary between stable configura- 
tions with N and N — 1 electrons). 

As mentioned in Section [IJ an external magnetic field 
can induce changes in the spin configuration even when 
the number of electrons confined in the QD remains con- 
stant. Slater's rule can be used only when there is a 
change in occupancy of one eigenlevel. If a change in 
the number of electrons is accompanied by spin rotations 
(that is, the occupancy of several levels is changed at the 
same time) , it becomes necessary to use several transition 
states sequentially (this increases the amount of compu- 
tations but does not invalidate Slater's rule). However, 
the number of such situations is small so that the over- 
all application of Slater's rule significantly decreases the 
amount of computational time required to calculate the 
full charging diagram in magnetic fields [27|. 

Having obtained /x(AT), we can determine the electron 
addition energy by first computing the so-called capaci- 
tive energy |23j ]: 

A 2 (N)=fx(N)-fx(N-l). (8) 

We evaluate this quantity at the bias Vg (N) correspond- 
ing to the addition of the N-th. electron to the system 
so that A2(iV) = — fi(N — 1). Together with the charg- 
ing voltage Vg(N), it is this quantity which is analyzed 
in detail in the following Section. Note also that this 
definition of the addition energy is consistent with other 
calculations where a rigid external potential is used to 
create electron confinement 0, 0, H, El [n| . 

IV. RESULTS AND DISCUSSION 

A. Addition energy spectra at zero magnetic field 

We first consider the addition energy spectra of the cir- 
cular and rectangular QDs in the absence of a magnetic 
field. The calculated spectra are shown in Fig. Q]together 
with the corresponding experimental spectra taken from 
Ref. 0. The spectrum of the QD in the circular mesa 
(Fig. ^ lower part) has been discussed in great detail 
earlier [fj so here we just briefly review it. It exhibits 
pronounced maxima for two and six electrons due to the 



first and second shell closures characteristic of QDs with 
parabolic 2D circular confinement. The peaks and val- 
leys are a consequence of the interplay between confine- 
ment and many-body effects. For two electrons in the 
QD, the lowest single-particle state E\ is fully occupied. 
The third electron populates the next available eigenstate 
with £2 > £i 5 thus making addition of the electron much 
more energetically costly than for the case of the second 
electron. The same situation is repeated for N = 6 when 
the second shell is complete. A smaller peak at N — 4 is 
due to the fulfillment of Hund's first rule. The total spin 
S of the iV-electron system is equal to zero (singlet) for 
N = 2 and 6 and to one (triplet) for N — 4. The agree- 
ment between the calculated and experimental spectra is 
very good, for both peaks and valleys. From the energy 
separation between the two lowest adjacent eigenvalues 
we can deduce that the confinement strength fku in the 
empty QD is about 6 meV. 

The analysis for the rectangular QD spectra is more 
challenging. From a simple model of a 2D anisotropic 
harmonic oscillator with eigenenergies £n x .n H = f^x(n x + 
1/2) + huj y (n y + 1/2), with quantum numbers n x , n y — 
0, 1, one can see that in general, except at certain spe- 
cial deformations, there are no degeneracies in the eigen- 
spectrum at T 6J. This means that the shell structure 
of a circular QD with strong maxima for 2, 6, 12, ... elec- 
trons in the addition energy spectra will not hold here, 
and one should expect maxima for most even electron- 
numbers. This simple interpretation is largely confirmed 
by our calculations. 

Experimentally, because of side wall etching, the sizes 
of the fabricated mesas can be determined only with 
about 5 % accuracy g|. To assess the effect of this un- 
certainty for the rectangular QDs, we first performed cal- 
culations of the addition spectra for three systems with 
geometrical aspect ratios equal and ±5 % of the value 
given for rectangular QD X: namely (3 — 1.306, 1.375 
and 1.444. Note that the last structure has an aspect 
ratio that is very close to the nominal value (3 = 1.44 
for rectangular QD Y [6J. Hence, we compare our cal- 
culated spectra to spectra for both rectangular QDs X 
and Y. The results are also given in Fig. ^ middle and 
upper part respectively. Surprisingly our calculations for 
f3 = 1.444 are in excellent agreement with the experimen- 
tal data for rectangular QD Y, whilst none of the 5 % 
variation-spectra are able to describe the experimental 
results as well for rectangular QD X, although the calcu- 
lated spectrum for the structure with j3 = 1.444 is clos- 
est to that measured. In fact, the experimental addition 
spectrum for X is similar in some respects to the one for 
the circular QD indicating that this structure may have 
an effective aspect ratio much lower than can be assumed 
on the basis of its geometrical aspect ratio alone. This 
conclusion is also corroborated by examining the mag- 
netic field dependence of the charging digrams for these 
structures (see Section HV _B|I . 

Further inspection of the calculated spectra reveals 
that the difference (contrast) between the peaks and val- 
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leys for the rectangular QDs is much lower than for the 
circular QD. We also note that the peak for N = 2 is 
significantly suppressed, and a 10 % difference in [3 gives 
rise to the N = 6 feature changing from a minima to a 
maxima (increasing (3 from 1.306 to 1.444). This shows 
that the electronic structure in rectangular QDs is very 
sensitive to even small distortions of the confinement po- 
tential and the precise balance between all contributions 
to the total energy is essential for its accurate descrip- 
tion. On the other hand, the four electron peak appears 
to be robust and is present in all the calculated addition 
energy spectra. This is because there is a trivial shell clo- 
sure at N = 4 since there is now no degeneracy between 
the second and third lowest single-particle eigenvalues as 
is the case for the circular QD. 

In order to better understand the electronic properties 
of a rectangular QD, we plot the wave- functions for the 
four lowest eigenvalues for the case of an empty QD at 
Vg — —1.6 and —1.3 V, and also for the QD contain- 
ing six electrons at Vg — —1.3 V (Fig. [2J). One can see 
that the " spread" of the wave- functions progressively in- 
creases from (a) to (c) not only because the presence of 
electrons in the QD induces a repulsive Coulomb poten- 
tial which effectively "flattens" the bottom of the conduc- 
tion band, but also this effect occurs in the empty QD by 
solely varying the gate bias. This is a consequence of the 
three-dimensionality of the mesa structure, which cannot 
be reproduced in model 2D calculations with a rigid con- 
finement potential. The doped regions just above and 
below the QD are never fully depleted, and changes in 
the gate bias affect the charge distribution in them as 
well. The resulting changes in the potential are carried 
over into the QD and lead to the effective decrease of 
the confinement strength as the gate bias becomes more 
positive. The same effect is also present in circular QDs 
It can be seen that the confinement decreases with 
increasing gate bias and the population in the QD. At 
the same time, the potential deformation ratio S also in- 
creases, i.e. the QD becomes more elongated. This is a 
manifestation of a self-consistent effect due to which elec- 
trons tend to localize as far as possible from each other 
in order to minimize the Coulomb repulsion. This results 
in a more effective decrease of the confinement along the 
x-direction than the y-dircction. 



B. Charging diagrams in a magnetic field 

In Fig. |31 we plot the calculated charging voltage for 
N = 1 to N = 5 electrons as a function of magnetic field 
(charging diagram) for the circular QD. Again these re- 
sults have already been extensively discussed elsewhere 
P-l24| but they are used here to compare the properties of 
circular and rectangular QDs. Comparison with exper- 
iment for the circular QD shows that overall agree- 
ment is very good, albeit the confinement potential at 
zero field in our model structure is somewhat weaker, i.e., 
the curves are more closely spaced. The charging volt- 



ages generally (in the considered range) increase with 
magnetic field since the effective confinement becomes 
stronger. The curves corresponding to the charging of 
N > 2 electrons exhibit "cusps" corresponding to vari- 
ous magnetic field induced spin and angular momentum 
transitions. For the N — 3 charging curve arising from 
the addition of the third electron to the two-electron sys- 
tem, the cusp near 5 T (depicted by a square) is due to 
a change in the ground state configuration of the two- 
electron system at B « 5.6 T (the singlet-triplet transi- 
tion) which affects the addition energy of the third elec- 
tron. The shift from 5.6 T to 5.0 T is due to a screening 
effect and is consistent with the experiment 7] . The cusp 
near 6 T (triangle) reflects an increase in the total spin 
of the three-electron system from S = 1/2 to S = 3/2, 
namely below this point two electrons are spin-up and 
one is spin-down while above this point all three elec- 
trons in the QD are spin-up and form a spin-polarized 
system. The dashed line in the 4 to 7 T region shows 
an excited state where three electrons have a total an- 
gular momentum L = 2 but total spin S = 1/2. In our 
calculations this state never becomes the ground state, 
which is different from recent results based on the config- 
uration interaction method pij ] where it was found that 
the transition from the ground state (L — 1, S = 1/2) 
to this state (L = 2, S = 1/2) was responsible for the 
lower field cusp. However, the predicted difference in the 
charging voltages is so small that given the large number 
of experimental parameters (such as dopant distribution, 
geometrical size of the structure, conduction band off- 
sets, etc.) whose values cannot possibly be determined 
precisely, the real cause for this feature is difficult to as- 
certain at best. 

In general, the rightmost cusp in the Fig. [^charging di- 
agram always corresponds to complete spin polarization 
of the electron system. The magnetic field at which the 
formation of this state occurs increases with the number 
of electrons since a stronger field is required to overcome 
the large kinetic energy accompanying single occupancy 
of consecutive orbitals. Cusps in the ./V-electron curves 
in the vicinity of B = 5 T are either due to changes in the 
(N — l)-electron configuration or due to single spin ro- 
tations in the ./V-electron system. The two cusps around 
B ks 0.25 T in the N — 4 and N = 5 curves mark the 
breakdown of Hund's first rule filling in the four-electron 
system and respective change in the addition energy of 
the fifth electron due to the decrease of exchange energy 
in the four electron system, similar to the N — 3 curve. 
The cusp at B w 1.5 T in the N = 5 curve can be un- 
derstood in terms of the Fock-Darwin spectrum describ- 
ing non-interacting particles Around that point, an 
electron undergoes a transition in the quantum numbers 
moving from the Fock-Darwin state with [n, I) = (0, —1) 
to (n, I) — (0, 2), thereby increasing the total angular mo- 
mentum of the system while keeping the total spin value 
constant. 

The calculated charging diagram of the rectangular 
QD with P = 1.444 shown in Fig. 0] superficially looks 
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similar to that of the circular QD. However, all cusps 
are less pronounced and the separation between consecu- 
tive curves is generally smaller. The former is due to the 
absence of a clear shell structure while the latter finds 
its origin in the decreased separation between the lower 
eigenlevels arising from the weaker confinement in the 
x-direction. We also observe that the magnetic fields at 
which transitions occur are approximately the same as 
in the circular QD (compare with Fig. [3J. This can be 
explained by noticing that the QD surface area and the 
resulting electron density in both QDs are approximately 
equal. Since in the DFT the electron density determines 
the magnitude of the Coulomb and exchange interactions, 
it takes roughly the same amount of energy in the two QD 
systems to change these interactions in order to induce 
single electron transitions. The charging curves move up- 
ward slightly faster with magnetic field than those in Fig. 
|3J The reason for this behavior is again, compared to the 
magnetic confinement lu c = eB/mc, a smaller value of 
hui x so that the magnetic field affects the related energy 
levels more strongly. Note that for the rectangular QD 
the N = 4 curve corresponds, at zero field, to the charg- 
ing of the fourth electron to form a singlet 5 = state 
(two electrons with spin-up, two with spin-down), and 
this behavior is clearly different from that for the circu- 
lar QD where four electrons obey Hund's first rule close 
to zero magnetic field. The singlet state continues to be 
the ground state up to B w 5 T where the S — 1 state 
(shown by a dashed curve) briefly takes over. 

At low magnetic fields (but above B w 0.25 T), the 
calculated charging voltages for both N = 3 and N = 4 
in the circular QD (Fig. |3J) clearly decrease whilst the 
corresponding curves for the rectangular QD are initially 
"flat" before moving upward (Fig. This is due to 

the different response of the second lowest single-particle 
state to the magnetic field. For a circular potential it first 
decreases and then increases with magnetic field, whilst 
a sufficiently large potential deformation results in the 
monotonic increase of this eigenvalue with magnetic field. 
This is what happens in our modelled structure where 
5 > 2 (see Fig. |5J). At small magnetic fields the experi- 
mental charging diagram for rectangular QD X [see Fig. 
EJa)] actually looks very similar to both the experimental 
and calculated diagrams for a circular QD (see Ref. and 
Fig. 13 respectively) , although notably the cusps related 
to Hund's first rule filling for N = 4 are absent near 
T. This indicates that the effective aspect ratio S char- 
acterizing the confinement potential in rectangular QD 
X is much closer to one than could be deduced from its 
geometrical aspect ratio j3. On the other hand, the mea- 
sured charging diagram for rectangular QD Y [Fig. Elb)] 
is closer in appearance to the calculated diagram for the 
rectangular QD shown in Fig. 01 I n particular, curves 
for N — 3 and N = 4 move gradually upwards, whilst 
the curve for TV = 5 moves only slightly downwards. The 
experimental charging spectra for both circular and rect- 
angular QDs nonetheless reveal a weaker dependence of 
the charging voltages on the magnetic field and a larger 



separation between consecutive curves (within 30%) not 
observed in the theoretical spectra |2^. Both these dis- 
crepancies may be explained if the confinement is indeed 
stronger (in both x and y directions) in the experimental 
devices, so that the magnetic field is unable to perturb 
it appreciably. 

C. Singlet-triplet transition in a two electron 

system 

One of the most striking differences between the cal- 
culated charging behavior of the QDs in the circular and 
rectangular mesa structures (Figs. |3]and0J) is found in 
the two electron system. The dashed curves for N = 2 
in both plots represent excited states. At low magnetic 
fields, below the cusp on the charging curve, the excited 
state is a triplet (S = 1) whilst at high magnetic fields, 
above the cusp, it is a singlet (S = 0). This cusp in fact 
marks the singlet-triplet (S-T) transition. From the full 
and dashed curves in Figs. |3and^ we can extract (see 
Eqs. 7 and[HJ the energy separation between the singlet 
and triplet states, namely 

J = J E T (2)- J B S (2) = AJ(2)-Af(2), (9) 

where symbols T and S stand for triplet and singlet states 
respectively. This singlet-triplet energy difference for a 
QD with two electrons is very interesting because of its 
potential importance for quantum computation. To the 
best of our knowledge this is the first calculation of this 
important characteristic in a realistic system taking into 
account the full 3D nature of the problem. Hitherto, 
the singlet-triplet separation has been considered in a 
realistic non-circular QD using 2D confinement potentials 
produced by the gate voltages only |25|. Our theoretical 
data for both circular and rectangular QDs as a function 
of magnetic field are shown in Fig. |SJ One can see that in 
the range of investigated magnetic fields the behavior of J 
in the two structures is radically different. In the limit of 
large magnetic field J — > for the rectangular mesa QD 
whilst it continues to become more negative without 
any saturation in the case of the circular QD. Also, the 
value of J for the rectangular QD is about three times 
smaller than for the circular QD at zero magnetic field. 

In order to gain a better insight into the calculated be- 
havior of J in a magnetic field, we show in Figs. {7\ and 
03 respectively the eigenfunctions of the occupied singlet 
and triplet states in a rectangular QD with j3 = 1.444. In 
Fig. 0we see that at low fields both of the electrons in the 
singlet state occupy the same space in the QD, i.e., their 
wave-functions are identical, and they become ever more 
localized at the center of the QD with increasing field. 
However, above B « 5.8 T the electrons start to localize 
in opposite regions of the QD along the x-direction induc- 
ing a symmetry "breaking" so that the wave- functions no 
longer have the spatial symmetry of the confinement po- 
tential around its minimum at the center of the QD. This 
effect appears similar to so-called Wigner crystallization, 
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i.e., spatial separation of electrons, caused by the applied 
magnetic field Q| . Further increase of the magnetic field 
gives rise to further localization of the electrons. In the 
case of the triplet state (Fig. |SJ, there is no such radi- 
cal reconstruction of the electron wave- functions. They 
just become more localized around the center of the QD. 
The symmetry breaking in the singlet state also results 
in a change of the total electron density. It attains two 
"humps" similar to the electron density in the triplet 
state. The reconstruction of the singlet electron density 
is responsible for the predicted evolution of the singlet- 
triplet splitting at large magnetic fields. When the elec- 
tron densities of both states become similar, exchange 
and Coulomb energies also have close values. Coupled 
with the fact that the single-particle energy separation 
decreases with increasing field (the energy levels condense 
in to the lowest Landau band Q), this results in a van- 
ishing J. Interestingly, as can be seen from Figs. [7| the 
overlap between the electron wave-functions at B = 8 
T is still appreciable (and close in value to that one at 
B = 6 T), and yet J — s- 0. 

The symmetry breaking does not occur in a circular 
QD where the electron density of the singlet state retains 
the symmetry of the confinement potential in the range of 
magnetic fields considered in this work. As a result, the 
exchange energy of the triplet state is always smaller than 
that of the singlet. Starting from zero field, the decrease 
in the separation between the two lowest single-particle 
energies as the magnetic field is increased leads to the 
triplet becoming the ground state at B « 5.6 T, but the 
singlet-triplet energy difference J does not exhibit any 
saturation behavior at higher field. 

The behavior of the calculated singlet-triplet separa- 
tion for the circular QD is in good agreement with ex- 
perimental data obtained by excited state spectroscopy 
which enables the observation of electron tunneling 
through both ground and excited states. Although there 
is extensive experimental data for the circular QD, de- 
tailed excited state spectra for the vertical rectangular 
QDs are currently unavailable. Clearly in light of the 
predictions concerning the behavior of J, it would be 
interesting to perform such experiments in the future. 
At present, we can only identify the two-electron S-T 
transtion for both rectangular QDs in the charging di- 
agrams of Fig. [S] which just show the evolution of the 
ground states. We do, however, note that recently, co- 
tunneling was utilized to extract the values of the singlet- 
triplet energy sepration for two electrons in a lateral QD 
with a spatially elongated (approximately harmonic) con- 
finement potential |2q ] where it was found that with in- 
creasing magnetic field, J changes sign and then satu- 
rates, consistent with our calculations. 



V. CONCLUSION 

In this work we studied the electronic properties of 
vertical quantum QDs in circular and rectangular mesas 
using a multiscale density-functional approach within the 
local spin density approximation. The QD was treated 
fully quantum mechanically and the surrounding regions 
were described by a semiclassical Thomas-Fermi approx- 
imation. The calculations were done on a 3D mesh by 
means of a finite element method. Slater's rule (the tran- 
sition state technique) was utilized to study charging of 
the QDs. With this approach, we were able to calculate 
the zero magnetic field electron addition energy spectra 
in the QDs as well as study the electron charging dia- 
grams in a magnetic field. 

We found that the zero field addition energy spectrum 
and charging diagram (charging voltages as a function 
of magnetic field) of a circular QD are in good agree- 
ment with available experimental results. For rectangu- 
lar QDs, we found that a 5 % variation from the nominal 
aspect ratio can give rise to marked differences in the 
addition energy spectra. On the basis of a comparison 
with experimental data in a magnetic field for rectangu- 
lar mesa structures (Rectangles X and Y in Ref. 0) with 
different aspect ratios, we conclude that rectangular QD 
X most likely has a smaller effective aspect ratio than 
previously reported, whilst our calculated results are in 
good agreement with those obtained for the rectangular 
QD Y. 

From the charging diagrams, we also extracted the cal- 
culated magnetic field dependence of the singlet-triplet 
energy separation J for a two electron system. Up 
to the singlet-triplet transition, this quantity decreases 
with magnetic field for both circular and rectangular 
QDs. However, beyond that transition, in the limit of 
a strong field, J for the rectangular QD asymptotically 
approaches zero whilst it continues to become more nega- 
tive in the case of the circular QD in the range of consid- 
ered magnetic fields. Analysis of electron wave-functions 
revealed that the behavior of the singlet-triplet energy 
separation in a rectangular QD can be attributed to the 
"breaking" of the symmetry of the singlet wave- functions. 



Acknowledgments 

This work has been supported by DARPA QUIST pro- 
gram DAAD19-01-1-0659 and Strategic Application Pro- 
gram (SAP) at NCSA. 



[1] K.K. Likharev, Proc. IEEE 87, 606 (1999). 
[2] D. Loss and D.P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 



[3] S. Tarucha, D.G. Austing, T. Honda, R.J. van der Hage, 
and L.P Kouwenhoven, Phys. Rev. Lett. 77, 3613 (1996). 
[4] S.M. Reimann and M Manninen, Rev. Mod. Phys. 74, 



8 



1283 (2002). 

[5] P. Matagne, J.-P. Leburton, D.G. Austing, and S. 

Tarucha, Phys. Rev. B 65, 085325 (2002). 
[6] D.G. Austing, S. Sasaki, S. Tarucha, S.M. Reimann, 

M.Koskinen, and M. Manninen, Phys. Rev. B 60, 11514 

(1999). 

[7] L.P. Kouwenhoven, T.H. Oosterkamp, M.W.S. Danoesas- 

tro, M. Eto, D.G. Austing, T. Honda, and S. Tarucha, 

Science 278, 1788 (1997). 
[8] A.V. Madhav and T. Chakraborty, Phys. Rev. B 49, 8163 

(1994); S. Tarucha, T. Honda, D.G. Austing, Y. Tokura, 

K. Muraki, T.H. Oosterkamp, J.W. Janssen, and L.P. 

Kouwenhoven, Physica E 3, 112 (1998); T. Ezaki, Y. 

Sugimoto, N. Mori, and C. Hamaguchi, Semicond. Sci. 

Technol. 13, Al (1998); K. Hirose, N.S. Wingreen, Phys. 

Rev. B 59, 4604 (1999); I.-H. Lee, Y.-H. Kim, and K.- 

H Ahn, J. Phys.: Condens. Matter 13, 1987 (2001); B. 

Szafran, F.M. Peeters, S. Bednarek, and J. Adamowski, 

Phys. Rev. B 69, 125344 (2004). 
[9] P.S. Drouvelis, P. Schmelcher, and F.K. Diakonos, Phys. 

Rev. B69, 155312 (2004). 
[10] Y. Tokura, S. Sasaki, D.G. Austing, and S. Tarucha, 

Physica B 298, 260 (2001). 
[11] P. Matagne, and J.-P. Leburton, Phys. Rev, B 65, 155311 

(2002). 

[12] P. Matagne and J.-P. Leburton, in Quantum Dots and 
Nanowires, Eds. S. Bandyopadhyay and H. S. Nalwa, 
ASP, California (2003). 

[13] J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 
(1992). 

[14] F. Ancilotto, D.G. Austing, M. Barranco, R. Mayol, K. 
Muraki, M. Pi, S. Sasaki, and S. Tarucha, Phys. Rev. B 
67, 205311 (2003). 

[15] S.M. Sze, Physics of Semiconductor Devices, 2nd Edi- 
tion, Wiley Interscience, 1981. 

[16] P. Matagne, and J.-P. Leburton, Phys. Rev, B 65, 235323 
(2002). 



[17] K.-J. Bathe, Finite Element Procedures in Engineering 
Analysis, Prentice-Hall, Inc., New Jersey, 1982. 

[18] S.R. White, J.W. Wilkins, and M.P. Teter, Phys. Rev. B 
39, 5819 (1989); J.E. Pask, B.M. Klein, C.-Y. Fong, and 
P. A. Sterne, ibid. 59, 12352 (1999). 

[19] S. Balay, K. Buschelman, V. Eijkhout, W.D. Gropp, D. 
Kaushik, M.G. Knepley, L.C. Mclnnes, B.F. Smith and 
H. Zhang, PETSc Users Manual, ANL-95/11 - Revision 
2.1.5, Argonne National Laboratory, 2004. 

[20] J. Fettig, D.V. Melnikov, N. Sobh, and J.-P. Leburton 
(unpublished) . 

[21] J.C. Slater, Adv. Quantum Chem. 6, 1 (1972). 

[22] L.R.C. Fonseca, J.L. Jimenez, J.-P. Leburton, and R.M. 
Martin, Phys. Rev. B57, 4017 (1998). 

[23] I.-H. Lee, K.-H Ahn, Y.-H. Kim, R.M. Martin, and J.-P. 
Leburton, Phys. Rev. B 60, 13720 (1999). 

[24] B. Szafran, S. Bednarek, and J. Adamowski, Phys. Rev. 
B 67, 115323 (2003). 

[25] J. Kyriakidis, M. Pioro-Ladriere, M. Ciorga, A.S. Sachra- 
jda, and P. Hawrylak, Phys. Rev. B 66, 035320 (2002). 

[26] D.M. Zumbiihl, CM. Marcus, M.P. Hanson, and A.C. 
Gossard, Phys. Rev. Lett. 93, 256801 (2004). 

[27] Our previous calculations flrjl have established the valid- 
ity of Slater's rule for the number of electrons N > 1. 
Nevertheless, for completeness we apply Slater's rule to 
determine the charging point for the first electron as well. 

[28] Surprisingly the smaller separation between the calcu- 
lated curves does not affect the addition energy spec- 
tra which are in good agreement with measured data 
(Section II V All . This is because the addition energy is 
defined not only by the gate voltage difference between 
consequitive charging points but rather by the slope of 
n{N, Vg) vs. Vg which in itself is a highly non-trivial 
function of the QD device parameters and electron num- 
ber (for discussion, see, e.g., [TM |24T| and references 
therein). 



°1 2 3 4 5 6 7 8 

Number of Electrons 

FIG. 1: Addition energy spectra for the circular QD (lower 
part), the rectangular QD X (middle part), and the rectangu- 
lar QD Y (upper part) at zero magnetic field. The parts are 
shifted by 3 meV in the vertical direction for clarity. In all 
parts solid lines represent results of calculations, while dashed 
lines stand for experimental data. Dotted and dash-dotted 
lines for rectangular QD X are the results of calculations for 
structures with f3 = 1.306 and 1.444 respectively. The cal- 
culated curve for the rectangular QD Y is the same as the 
dash-dotted curve for the rectangular QD X. 
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FIG. 2: Two-dimensional contour plots of the wave-functions 
for the four lowest eigenvalues (from top to bottom) in the 
rectangular QD with (3 = I AAA: (a) N = 0, V G = -1.6 V, (b) 
N = 0,V G = ~1.3 V, (c) N = 6, V G = -1.3 V. Coordinate z 
is fixed in the center of the QD. 
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FIG. 3: Calculated charging diagram (solid curves) for ground 
states in a magnetic field for the circular QD. The dashed lines 
show some excited states. Triangles (?) define boundaries of 
regions with different total spin and angular momentum in 
the ground state of the A^-electron system. Squares (■) show 
features in the ground state of the iV-electron system due to 
reconstruction of the (N — l)-electron system. 
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FIG. 4: Same as in Fig. EDbut for the rectangular QD with j3 = 
1.444. In this case triangles (?) define boundaries of regions 
with different total spin and (n x , n y ) quantum numbers in 
the iV-electron system. 
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FIG. 5: Experimental charging diagrams showing ground 
state evolution with magnetic field for (a) rectangular QD 
X, and (b) rectangular QD Y. The position of the N = 2 
singlet-triplet (S-T) transition is marked. 




FIG. 6: Calculated singlet-triplet energy separation J in the 
two-electron system as a function of the magnetic field. /3 = 
1.444 for the rectangular QD. 
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FIG. 7: Single-particle electron densities in the x-direction 
in the rectangular QD with f3 = 1.444 for the singlet state 
at different magnetic fields (coordinates y and z are fixed in 
the middle of the QD). Arrows (j) and (J.) correspond to 
the electrons with spin up and down. The data for B = 
0, 4 T and 6, 8 T are obtained at V G = -1.45 and -1.47 
V respectively. The curves for different magnetic fields are 
shifted in the vertical direction by 2.5 x 10 23 cm -3 for clarity. 
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FIG. 8: Same as in Fig. Qbut for the triplet state. Curves 
with circles (•) and pluses (+) are for the ground and first 
excited occupied states respectively. The curves for different 
magnetic fields are shifted in the vertical direction by 2.0 x 
10 23 cm -3 for clarity. 



